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43 Abstract 

A phenomenological method is developed to consider elastic guided wave 
propagation in complex curved waveguides. The theory on guided wave prop- 
agation in hollow circular cylinders is used in order to verify the method. The 
results are compared with solutions obtained by the Helmholtz decomposition 
method. The method given here is used to derive new analytical relations 
for torsional and longitudinal phase velocity dispersion curves. In addition, 
a new excitation method is proposed and examined for flexural modes. Some 
conceptual results and applications of the method are discussed. 
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1. Introduction 

Ultrasonic guided waves are utilized for rapid, long-range inspections of 
various structures even when access to the structure is limited Guided 
waves are used extensively for inspecting isotropic single layer plates 0] , com- 
posite plates [ij], pipelines 0,0], rail roads [f|, adhesive joints Q, and more. 
In addition, current industrial demand for applications of guided waves is 
growing. Often, new applications of guided waves involve complex struc- 
tures. 

Guided waves are more complicated than bulk waves because of their disper- 
sive nature, infinite number of propagating modes, wave velocity and wave 
structure variations with mode and frequency, and waveguide boundary con- 
ditions SB [Tof . Despite the complex nature of guided waves, a complete 
understanding of guided wave propagation in a waveguide is necessary for 
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robust theoretically driven inspection. For example, a prior knowledge of 
group velocity, wave structure, and skew angle of excited modes is required 
to estimate the location and type of a defect in a composite plate. On the 
other hand, if the complexities of guided waves are quantitatively understood 
and utilized, the sensitivity and penetration power of guided waves can be 
significantly improved for inspection purposes . For instance, wave propaga- 
tion characteristics of non-axisymmetric flexural modes in a pipe are used for 



synthetic and active focusing [111 . 1 121 ] . Another example of how theoritically 



understanding of guided waves can be applied to improve inspections is using 



wave structure matching for inspection of adhesive joints HI • 
Guided wave dispersion curves are crucial for presenting the phase and group 
velocity variation with frequency. Furthermore, dispersion curves reveal the 
dispersivity characteristics of an excited mode. This information is useful for 
signal post-processing. Dispersion curves are also used to determine suitable 
modes for excitation, the angles of angled beams, and the spacing of comb 
type transducers to excite special modes 

While dispersion curves, wave structure, mode selection, scattering proper- 
ties, and wave propagation characteristics of guided waves should be stud- 
ied for efficient and successful inspection, these concepts can be considered 
in a straight forward fashion for some simple geometries like single layer 
plates, multilayer plates, and pipes, where Navier's equation can be solved 
by Helmholtz decomposition or by the partial wave method [To[. Finite 
element analysis can be used when considering complex waveguides, but the 
method suffers from a fundamental inability to provide a generic solution 
for further sophisticated implementation of guided waves. In addition, com- 
putational cost and time consumption of finite element analysis is high, so 
it cannot be used for geometrically large waveguides. Semi-analytical finite 
element analysis (SAFE) can decrease the computation cost but is applicable 
for structures only with constant, albeit arbitrary cross section 14 . 
On the other hand, understanding wave propagation properties in some im- 
portant but complex curved geometries is highly desirable. For example, 
elbows are an essential element of pipelines, so an effective long range in- 
spection of a pipeline depends strongly on the ability to inspect beyond 



elbows [15]. The understanding of wave propagation in an elbow should al- 
low data interpretation of reflected waves from beyond the elbow. A growing 
number of NDE projects are becoming of interest in which complex curved 
surfaces are involved. Analytical methods like Helmholtz decomposition and 
the partial wave method cannot provide solutions for geometrically complex 
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curved structures; and applicability of numerical solutions is also limited. 
Therefore, developing a practical method to consider guided wave propaga- 
tion in complex curved waveguides would be highly beneficial for extending 
the applications of guided waves for inspection purposes. 
In this work, a physical understanding of guided waves propagation is pre- 
sented. A method for predicting guided wave propagation in complex waveg- 
uides is developed by considering the current understanding of different as- 
pects of wave propagation. Then, upon the new understanding given here, a 
quantitative method is developed. The method is used to derive an analytical 
relation for torsional and longitudinal phase velocity dispersion curves in a 
hollow circular cylinder. The results of the method are compared with other 
available solutions. Limitations of the method are considered in section HJ 
In section O a new method for flexural mode excitation is proposed based 
on a newly developed perspective of guided wave propagation. Finite Ele- 
ment models are used to show that the method effectively excites the desired 
flexural modes. A summary of the method and potential applications are 
considered later in the text. 



2. Method Description 

A method is developed to study wave propagation in complex curved 
waveguides. There are a few approximation methods that can simplify the 
wave equation by low or high frequency approximations. They are based 
upon a Taylor series expansion of the wave equation and then eliminating 



the higher orders based on some assumptions [16|, |17|, [18|, |19| . In this work, a 
phenomenological approach is adopted instead of a mathematical approach. 
It is shown that the method is accurate and simple to apply in practical 
situations. The method is motivated by the following facts. 

1. A plane wave in an infinite medium is described by a specific direction 
of propagation and also a sinusoidal change in amplitude which is a 
function of time and distance along the wave propagation direction. 
One wavevector is chosen on behalf of an infinite number of wavevec- 
tors that exists all around the propagation region, because all of the 
wavevectors have the same direction and amplitude in the case of plane 
wave. Consider the generalized case where the wave propagation di- 
rection is different at each point. Then the plane wave symmetry is 
broken and the wave should be determined by a variable vector field 
that is a natural generalization of constant wavevectors. 
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2. The plane wave in infinite medium moves on a straight line, because 
the medium is assumed to be isotropic and flat. Consider a plane wave 
that moves in a two dimensional flat plate. The propagation path of 
the plane wave is considered as a straight line, but if the plate is curved 
the propagation path changes to a curved line. 

Therefore, in order to extend a plane wave propagation in an isotropic 
medium to a multi-directional wave propagation on a two dimensional curved 
surface, the above facts lead us to the following perspective. 

1. An arbitrary wave is defined by a vector field that is changing point 
by point according to the underlying geometry and initial emission 
conditions. 

2. An infinite number of paths can be defined on the surface that are 
tangential to the propagation direction vector field at each point. These 
paths are called ray-paths. In fact, ray-paths on a surface are geodesies 
of the surface. Geodesies are defined as straight lines on a curved 
surface. The geodesies of a surface can be determined by the metric of 



the surface and the following relation [29 

dx x > dx^ 1 dx 



+ r A = o (i) 

ds 2 + ^ ds ds ' 1 ' 

where s is affine parameter, A, [A, v — 1,2 and T^'s are the Christoffel 
symbols of the metric. 



3. The intensity of the wave will change sinusoidally along each ray-path. 

We call the method a ray-method. The method can be summarized as fol- 
lows. Assume that the excitation from a certain source on a complex curved 
surface is studied. Infinite emitted rays are considered where the excitation 
conditions determines the initial position, initial direction, and maximum 
intensity of each ray. Then, the geodesic equation determines the propa- 
gation path of each ray on the surface. The intensity of each ray changes 
sinusoidally along the ray-path. In order to find the intensity of the wave at 
each point, a superposition of all ray-path intensities that pass through the 
point at the same time is considered. 

In order to develop a complete method of studying guided wave propagation 
in a curved waveguide, though, the effect of thickness on guided waves should 
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be addressed. Until now, we have used rays to consider guided wave propa- 
gation. In particular, rays that lie on the surface of the waveguide. In this 
model, the thickness of the waveguide cannot influence ray propagation. In- 
deed, the ray-method cannot capture guided waves properties caused be the 
thickness of the waveguide like dispersion curves and wave structures. Now, 
the dispersion curves of a plate are compared with the dispersion curves of 
a pipe in order to get some insight on how the thickness can be included in 
the ray-method. 

Two phase velocity dispersion diagrams are defined for a plate: shear-horizontal 
dispersion curves (SHDCs) and Rayleigh-Lamb dispersion curves (RLDCs). 
SHDCs are derived from the following analytical relation [§[ 



= -7===F- ( 2 ) 



2or(fd) 



y/4(fd) 2 -v*<* 



T 
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Where ct, fd, and c p are the shear wave velocity, the frequency-thickness 
product, and the phase velocity of the nth mode, respectively. RLDCs are 
derived numerically from the Rayleigh-Lamb transcendental equation [gj]. 
Governing transcendental equations for deriving dispersion curves in a pipe 



are more complicated than for a plate [20|. The complicated transcenden- 
tal equations need to be solved in order to find torsional dispersion curves 
(TDCs) and longitudinal dispersion curves (LDCs) in a pipe. Figure Q] shows 
the axisymmetric part of TDCs and LDCs in a steel pipe (red) with s = 0.25 
and SHDCs and RLDCs in a steel plate (blue) on the same graph. Con- 
sidering the dispersion curves in Figure [T] reveals that the axisymmetric part 
of pipe dispersion curves nearly are the same as plate dispersion curves for 
higher frequencies. The primary differences are in the low frequency limit 
where L(0, 1) and L(0, 2) experience a jump to higher phase velocities. This 
jump for the pipe with s = 0.25 occurs around 25[kHz}. The frequency 
where this jump occurs tends toward zero as the s factor decreases. In most 
practical pipes the s parameter is less than 0.25, so we can approximate the 
axisymmetric part of pipe dispersion curves by plate dispersion curves. 
An axisymmetric mode in the context of the ray-method is represented by 
rays that move along the pipe axial direction. Thus, the similarity between 
the axisymmetric part of the pipe dispersion curves and the plate dispersion 
curves means that we should get plate type dispersion in a pipe, if we con- 



s is the ratio of pipe thickness to pipe mean radius 
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Plate Phase Velocity DC vs. Pipe Phase Velocity DC (s=0.25) 




50 100 150 200 250 300 350 
Frequency[kHz] 



Figure 1: Figure shows the axisymmetric dispersion curves in a steel pipe with s = 0.25 
(red) and the dispersion curves in a steel plate (blue). The diagram shows that the 
axisymmetric dispersion curves in the pipe are nearly the same as the dispersion curves 
in the plate for higher frequencies. 

sider a set of rays moving along the pipe axial direction. This observation 
leads to the notion that instead of using a single ray confined to the surface 
of a waveguide we could consider a curved plate spanning the thickness of the 
waveguide where the plate falls along a ray-path and is normal to both free 
boundaries of the waveguide. In this context, the axisymmetric dispersion 
curves of a pipe will be same as the dispersion curves of a plate. Hereafter, 
this concept will be referred to as ray-plate method. 

In the following section, guided wave propagation in a pipe is studied by the 
ray-plate method. New analytical relations are derived for pipe dispersion 
curves using the ray-plate method. Results from the ray-plate method are 
compared with exact results derived by the global matrix method. It is shown 
that the ray-plate method predicts accurate results for the pipe dispersion 
curves. In addition, it provides new insight on pipe flexural modes. 
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3. Verification of Ray-Plate Method 



In this section, the ray-plate method is applied to hollow circular cylin- 
ders. Guided waves in hollow circular cylinders have been studied extensively 



20[. The global matrix method with numerical calculation are often used 
in a complicated form to derive the pipe dispersion curves There are no 
analytical solutions for the torsional and the longitudinal dispersion curves 
of pipe. The ray-plate method is presented in detail and its results are com- 
pared with existing solutions for guided waves in hollow circular cylinders. 
As mentioned before, comparing dispersion curves between plates and pipes 
show that SHDCs in plates are approximately the same as the axisymmetric 
TDCs in pipes and also that the RLDCs in plates approximate axisymmetric 
LDCs in pipes. Also, the axisymmetric modes in the context of the ray- 
method corresponds to an infinite set of rays that move parallel to the pipe 
axial direction. This leads to the extension of the ray-method from mem- 
branes to thick-walled shells, if instead of a set of rays, a set of plates across 
the thickness of the surface is considered. Each plate propagates along a ray- 
path and is perpendicular to both boundary surfaces of the shell. Each plate 
carries Lamb waves and/or SH waves. These plates are called ray-plates. 
The ray-plate method claims that axisymmetric modes are ray-plates that 
move parallel to the axial direction of pipe. Ray-plates that move at a rel- 
ative angle with the axial direction can naturally become a candidate for 
explaining flexural modes. Let's construct dispersion curves in pipe by the 
ray-plate method base on the above insight about axisymmetric and flexural 
modes. 

As mentioned above a flexural mode is constructed by a set of ray-plates 
carrying Lamb waves or SH waves and moving in a direction that makes 
an angle a with the pipe's axial direction (z-direction). Consequently, the 
phase velocity of the flexural modes and corresponding axisymmetric modes 
are the same at a fixed angular frequency u. The only difference is that these 
modes are moving in a different direction than axisymmetric modes, namely 
by the angle a. Dispersion curves in a pipe are defined by choosing the z-axis 
as a preferred direction. The phase velocity of each mode is determined by 
c p = u)/k z , where k z is the z-component of the wavevector. Therefore, the 
phase velocity of a flexural mode can be derived by the following relation: 

M) = ^ = 1 >) (3) 

p kcos{a) cos(a) p ' 1 ] 
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where c p a ^ is the phase velocity of the axisymmetric mode at each frequency 
lo and a is the angle of wave propagation . The cylindrical shape of a pipe 
imposes a periodic displacement continuity boundary condition in the cir- 
cumferential direction. A flexural mode is defined as a Lamb or SH wave 
existing in a ray-plate, where the ray-plate makes angle a with the z-axis. 
So, the displacement field components can be written as 

u(r,<j),z) = fjM e i i kcos ( a ) z + k8ia ( a ) Iie f>- UJt ) , (4) 

In a pipe, if <fi is changed to <fi + 2ir, the displacement vector should be 
the same, meaning that u(r, <ft, z) = u(r, <ft + 2tc, z). This periodic condition 
imposes that 

ksm(a)R = m m — 0, ±1, ±2, (5) 

where R, k, and to are the mean radius of the pipe, the wavenumber, and 
the flexural order. From fl3]) and fl5]), the following relation is established 
between the phase velocity of a flexural mode and the phase velocity of its 
corresponding axisymmetric mode at a given frequency, 

C { v f) = ; 1 C j a l (6) 

p V^iW p 

The relation can be rewritten as 



(0,n) 

c^\fd) = Cp (7) 



(0,n) 
msCp 



_ 2n(fd) 

where Cp" 1 '^ is phase velocity of a mode of family n and flexural order to, fd is 
the frequency-thickness product and s = d/Ris the ratio of the pipe thickness 
to the pipe mean radius. Relation (I7j) imposes the following condition on the 
maximum number of flexural modes that can exist at each frequency (f), 
Lamb mode phase velocity c p °' n ^ , and radius of the pipe (R): 

m max = mt I -^y I (8) 

Eq.flE]) gives a new limitation on the number of flexural orders at each fre- 
quency. This relation can not be extracted simply from the global matrix 
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method and needs further consideration to be recognized. 
Eq.(J7]) is used to derive LDCs in a pipe. Mean radius and thickness of pipe, 
plate phase velocities at each frequency, and flexural orders should be sub- 
stituted in eq.© in order to drive the LDCs. Figure 2 shows a comparison 
between the LDC derived by the ray-plate method and the global matrix 
method. It is observed that at lower frequencies, flexural modes are sepa- 



rated more from each other than at higher frequencies. Figures 2(a) and 2(b) 
show that the derived LDCs in a 3. 5 [in] schedule 40 pipe by the ray-plate 
method are in good agreement with dispersion curves derived by the global 
matrix method. The following analytical relation for TDCs, including flex- 
ural modes, can be derived by using equation ([7]) and the shear horizontal 
dispersion relation (J2J): 

>*0(/d) = — 2*C T (fd) 

P y/4TT 2 {fd) 2 - (vr% 2 + m 2 s 2 )C T 

where Ct, fd, n, and m are shear wave velocity, frequency-thickness product, 
family order, and flexural order, respectively. Figure [3] shows a comparison 
between dispersion curves of a 3 [in] schedule 40 steel pipe derived by the 
ray-plate method and the global matrix method. The figure shows that the 
analytical relation (j5J) , derived by the ray-plate method, provides an accurate 
approximation for TDCs in pipes. 

It is simple to check that if the limit R — >• oo is considered, we get the shear 
horizontal dispersion curve relation. In addition, considering limit (fd) — > oo 
gives Cp 1 ' 11 ^ = Ct- Therefore, equation ([9]) is a consistent equation in those 
limits. In addition, equation fl2]) imposes the following relation to existence 
of torsional mode (n, m) as 



«V + 4m 2 <^- (10) 

The ray-plate method gives some insight on why flexural mode group veloc- 
ities are less than the axisymmetric mode group velocities. In the ray-plate 
method, flexural modes are related to rays whose propagation direction make 
an angle a with the z-axis. While the group velocities for all rays are the 
same, the propagation paths for the flexural modes will be longer than those 
of the axisymmetric modes. Therefore, the measured group velocities for 
flexural modes along the z-axis are less than those of their corresponding ax- 
isymmetric modes. Higher order flexural modes correspond to longer paths 
and hence lower group velocities. 



9 



Longitudinal Phase Velocity Dispersion Curve in a 3.5 [in] schedule 40 steel pipe 
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(a) 

Longitudinal Phase Velocity Dispersion Curve in a 3.5 [in] schedule 40 steel pipe 
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Figure 2: LDCs derived by ray-plate method and by global matrix are compared in these 
diagrams, a) compares the result of the ray-plate method and the global matrix method 
for 3.5[in] schedule 40 steel pipe , b) shows a magnified view of fig. (a) for lower frequencies. 
The figures present an excellent agreement between results of the ray-plate method and 
the global matrix method. 
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Torsional Flexural Phase Velocity DC in a 3 [in] sch. 40 steel pipe 




Frequency [kHz] 

Figure 3: This figure presents a comparison between the results of the ray-plate method 
and the global matrix method. It demonstrates that the ray-plate method closely approx- 
imates torsional dispersion curves for 3 [in] schedule 40 steel pipe. 

4. Ray-Plate Method Limitations 

Limitations of the ray-plate method for the prediction of guided wave 
propagation in complex curved surfaces are now considered. It is shown that 
the ray-plate method presents a relatively accurate approximation when the 
thickness to curvature radius ratio of a structure is less than 0.062 at each 
point of the structure. Otherwise, some deviations from the exact solutions 
will be observed. In many practical situations, however, the ratio of thickness 
to curvature radius is less than 0.062. Therefore, wave propagation in many 
practical curved surfaces can be considered by the ray-plate method. 
Considering the derivation of pipe dispersion curves by the ray-plate method 
in the previous section showed that in the context of ray-plate method dis- 
persion curves in a pipe can be derived from the dispersion curves of a plate. 
The ray-plates in general can be curved, for example, while axisymmetric 
modes in a pipe can only be constructed by flat ray-plates, the construction 
of flexural modes are possible only by using curved ray-plates. The curvature 
radius for the curved ray-plates in a pipe decreases for higher order flexural 
modes. The curvature radius for the axisymmetric ray-plate is infinite and 
for the flexural mode m = m max is the radius of the cylinder. 
For the sake of simplicity, the ray-plate method was based on flat ray-plates. 
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Therefore, any deviation from a flat plate because of waveguide geometry 
can potentially affect the predictions of the ray-plate method. Here, phase 
velocity deviations of SHDCs and RLDCs in a curved plate from the flat 
plate dispersion curves are considered in detail. Then, we calculate an ap- 
proximate turning point in which curvature of a curved ray-plate causes a 
significant deviation from the flat ray-plate and thereby invalidates the ray- 
plate method. A similar study has been done by Shao et al. 2l| which 
is limited to the consideration of SH waves. Other aspects of the problem 



have also been studied by Viktorov [22| and Qu et al. (23| under the topic 
of guided circumferential waves in a circular annulus. They showed that the 
dispersion curves in a curved plate experience some deviations from flat plate 
dispersion curves; and that the deviations are higher for a smaller curvature 
radius. 

A cylindrically curved surface (cylinder) is assumed to be long enough as 
to neglet edge effects. Guided waves are assumed to propagate in the cir- 
cumferential direction. Navier's equation is solved in cylindrical coordinates. 
Then, SH and Rayleigh-Lamb dispersion curves are derived for a variety of 
thickness to curvature radius ratios. 



Figures 4(a) and 4(b) show SH dispersion curves and Rayleigh-Lamb dis- 
persion curves for plates with different thickness to curvature radius ratios. 
These figures reveal that for value s < 0.062, the dispersion curves can be 
approximated by the plate dispersion curves. The s value for most of stan- 
dard pipe sizes is less than 0.062. Therefore the ray-plate method can be 
applied accurately for most of them. The situation is the same for many 
other structures that do not have small curvature radii. 



5. Flexural Mode Excitation 



Flexural modes in pipes have been known about for decades [20|, [24J, |26 
but only basic physical understanding of these waves can be extracted from 
the direct solution of Navier's equation by the Helmholtz decomposition 
method. The ray-plate method provides a new physical understanding of 
these waves. The method reveals that flexural modes are basically like ax- 
isymmetric modes which an angle relative to the axial direction 

in a pipe. This new understanding of flexural modes suggests a new method 
of exciting pure flexural modes in pipe. In this section, the new flexural mode 
excitation method is proposed and examined by finite element analysis. 



Flexural mode excitation has been considered by Li et al. 24] and Sun et 
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SH Phase Velocity Dispersion Curves 
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Figure 4: Figures show dispersion curves for plates with different ratio of thickness to 
curvature radius (s). a) shows SHDCs when s varies between 0.25 (cyan) and s p i a te = 
(black), b) shows RLDCs when s varies between 0.125 (cyan) and s p i ate — ( green 
square). The diagrams show that for s < 0.062, the dispersion curves of a curved plate 
can be approximated by the flat plate dispersion curves. 
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al. 26], where partial normal loading and angle beam wedges are used to 
excite flexural modes. Those works used normal mode expansion to find 
the amplitude of the excited axisymmetric and flexural modes with respect 
to the characteristics of the load. While flexural modes can be excited by 
these methods, they suffer two disadvantages. First, one pure flexural mode 
cannot be excited by these methods, several other flexural modes and also 
axisymmetric modes are excited at the same time. Next, the excited flexural 
modes can only partially cover the surface of a pipe, so a complete scan of 
the pipe is difficult with these waves. 

The active focusing technique is utilized for focusing a significant portion of 
guided waves energy at a desired point on the circumference position of a 



pipe. Li et al. [25] and Sun et al. [26|, |27|, |28j used appropriate superposition 



of angular profiles of flexural modes at each axial distance in order to focus 
at the desired circumferential angle. Synthetic focusing is also used for de- 
fect imaging in pipe. Quasi-axisymmetric waves are excited in the pipe and 
then the reflected flexural and axisymmetric modes are analysed to estimate 



the location and size of a defect. Hayashi et al. [TJJ developed a method 
to analyse reflected axisymmetric and flexural modes. They reconstruct an 
image of the defect by a time-reversal technique. Therefore, analyzing and 
implementing flexural modes is a crucial element of active and synthetic fo- 
cusing methods of guided waves in pipe. 

It has been shown that the simultaneous excitation of a ring of transducers 
around a pipe can excite quasi-axisymmetric modes very effectively. As we 
mentioned before, in the ray-plate method, flexural modes are the same as 
axsisymmetric modes, but they move at an angle corresponding to their flex- 
ural order. Considering this concept of flexural mode propagation suggests 
an angular excitation of the pipe can excite a pure flexural mode. Therefore, 
a normal load helical excitation (figJS]) is proposed to excite a special flexural 
mode. The helix angle is determined by the ray-plate method. Different 
angles of the helix can excite different flexural orders. The helix angle a is 
determined by using eq. (J5]) as follows: 

<«= <"> 

where / is frequency of excitation, Cp is the axisymmetric phase velocity of 
the desired family, and m is the order of the flexural mode that is excited. 
Flexural mode excitation in a 3.5 [in] schedule 40 steel pipe is considered. 
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Table 1: L(m, 1) and L(m, 2) propagation directions predicted by the ray-plate method 



f =50 [kHz] 


m=0 


m=l 
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6.27 


12.61 


19.12 


26.0 


33.1 
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24 


54.38 









Table 2: L(m, 1) and L(m, 2) phase velocities predicted by the ray-plate method 



f =50 [kHz] 


m=0 


m=l 


m=2 


m=3 


m=4 


m=5 


<(m,l)[ m /s] 


1525 


1534 


1563 


1614 


1695 


1820 


^L(m,2)[ m / S } 


5676 


6212 


9746 









Where the external diameter of the pipe is D = 3.5[m] and its thickness is 
d = 0.22[m]. Because of a clear separation of the flexural modes at lower 
frequencies, the frequency of excitation is chosen to be as low as 50[kHz]. 
Considering the longitudinal dispersion curves of the pipe shows that at this 
frequency two axisymmetric modes L(0, 1) and L(0, 2) exist where their phase 
velocities = 1 525 [m/s] and c p = 5676[m/s], respectively. Substitut- 

ing these data in eq. ffTTT) gives table CD of the ray-plate method predictions 
for the angle of propagation/excitation for m = 0,1, 2, 3, 4, 5. Furthermore, 
phase velocities are predicted as table [2] by the ray-plate method. Finite 
element analysis is utilized for considering the flexural mode excitation by 
the proposed method. A helically normal edge load is applied on a 3.5 [in] 
schedule 40 steel pipe (see figJSJ). The normal to the helix makes an angle 
of 20[cfeg] with respect to the z-axis. Table [1] shows that two flexural modes 
L(3, 1) and L(l,2) have propagation angle a = 19.12[deg] and a = 2A[deg], 
respectively. These angles are close to the excitation angle 20[cfeg]. Therefore 
we expect to observe both L(3, 1) and L(l, 2), but higher intensity at L(3, 1) 
with phase velocity 1624 [m/s]. A 50 [k H z] Gaussian modulated signal is used 
as an excitation signal. 

Figures 6(a) and 6(b) show excited waves in the pipe at t — 0.26 [msec] 
following the initiation of the excitation. The propagation direction of the 
wave clearly makes an angle with the z-axis. In addition, the wave is di- 
vided into two parts because of a difference in the group velocities. Figure 



6(b) shows the unwrapped version of the excited waves in the pipe. Phase 
velocities and flexural orders of these waves can be derived by taking the 
Fourier transform of the wave displacement field components on the z-axis 
and the pipe circumference. The z-component of the displacement field on 
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Figure 5: The new method of flexural mode excitation is presented in this figure. A helical 
edge load is applied on the pipe. The helix makes angle 20[deg] with z-axis. So we expect 
to have a relatively pure excitation of flexural mode £(3, 1). 



the external surface of the pipe is recorded. Taking the Fourier transform 
of these data on the z-direction gives us the wavenumber k. Then, using 
c p = u/k, the phase velocity of the wave around the pipe can be calculated. 
Figure 7(a) shows that there are two distinct wave numbers k = 193. 7[m -1 ] 
and k = 49.2[m~ 1 ] corresponding to c p = 1621[m/s] and c p = 6385[m/s], 
respectively. In the next step, the Fourier transform of the data is taken 



in the circumferential direction to extract the flexural orders. Figure 7(b) 
shows the order of the excited flexural modes at all axial distances. Figure 8 
shows the order of the excited flexural modes at z = 0.65[m] and z = l[m]. 
Figure 8(a) shows that at z — 0.65[m] the flexural mode is of order three. 

shows that at z — l\m] a flexural mode of order one has been 



Figure 
excited 



5(b) 



Considering the derived phase velocities for the excited modes in fig. 



7(a) 



the presented flexural order in figs. 7(b) , 8(a) , 8(b) along with tables [T] 
and [2] suggest that the flexural modes L(3, 1) and L(l,2) have been excited 
effectively. In conclusion, the FEM proves that the proposed method for 
flexural mode excitation works quite well. This implies that the interpreta- 
tion of flexural modes by the ray-plate method is correct. In addition, the 
predicted angles for flexural mode propagation by the ray-plate method are 
approximately correct. Therefore, the results in this section demonstrate the 
correctness and usefulness of the ray-plate method once more. 



16 



t=2.6( 


:-4[sec], Load=20[deg], 3.5[in] sch40 (m) 
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(a) 



Excited Guidedwave in 3.5[in] sch40, (20[deg] excitation angle) 




0.5 0.6 0.7 0.8 0.9 1 1.1 1.2 1.3 

z-axis distance [m] 



(b) 

Figure 6: Figures show the excited wave in a 3. 5 [in] schedule 40 steel pipe by 20[deg] helical 
excitation, a) shows directional propagation of the wave in the pipe, b) shows unwrapped 
view of the propagating wave. 
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Wavenumber at z-direction (k), 3.5[in] sch40,a=20[degree] 




Circuferebtial Angle [degree] 




Figure 7: Figure (a) shows the wavenumber k of excited modes. There are two distinct 
wavenumbers k = 193.7[to _1 ] and k = 49.2[m _1 ] corresponding to c p — 1621[m/s] and 
c p — 6385[m/s], respectively. According to tableQ, these phase velocities correspond to 
L(l,3) and L(2,l). 
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x Flexural Mode Content,(3.5[in]Sched40, a=20[degree]) 
61 1 1 1 r 




1 2 3 4 5 



Flexural Order 

(a) 

x 1 q-7 Flexural Mode Content,(3.5[in]Sched40, a=20[degree]) 
9i 1 1 1 r 




Flexural Order 



(b) 

Figure 8: Figure shows the order of excited flexural modes at z = 0.65[to] and z = l[m] 
at t = 0.26 [msec], a) shows that at z — 0.65[m] flexural mode has order three, b) shows 
that at z = l[m] a flexural mode of order one has been excited. 
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6. Summary 

A new method was developed for the consideration of wave propagation on 
complex waveguide surfaces. The method was derived from a phenomenolog- 
ical point of view and is computationally simple enough to apply to complex 
but practical situations. The method was extended to thick-walled waveg- 
uides where the thickness of the structure partially determines the dispersion 
and propagation properties of the guided waves. Pipe dispersion curves were 
used to prove the accuracy and the applicability of the method. For the first 
time, approximate analytical relations were derived for torsional and longitu- 
dinal dispersion curves. A physical interpretation of flexural modes in pipes 
was presented. Based on this new understanding of flexural modes, a helical 
excitation method was proposed for flexural modes. Finite element anal- 
ysis was utilized to show that helical excitation works appropriately. FEM 
showed that as the ray-plate method predicts, a 20[cfeg] helical excitation can 
effectively excite L(3, 1) and L(l, 2). In addition, limitations of the ray-plate 
method were considered for structures with the small radius of curvature. It 
was shown that the predictions of the method are accurate if the ratio of the 
thickness to the mean radius is less than 0.062. 

Continue this work, the proposed flexural mode excitation method is under 
study from several aspects. First, alternative methods to helical excitation 
are considered like using time delayed phased array. Also, the ability of ray- 
plate method to excite pure flexural modes at a wide range of frequencies 
and phase velocities is considered. 
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